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Abstract I discuss two of the possible sources of bi- 
ases in the determination of the IMF: binning and the 
existence of unresolved components. The first source 
is important for clusters with a small number of stars 
detected in a given mass bin while the second one is 
relevant for all clusters located beyond the immediate 
solar neighborhood. For both cases I will present results 
of numerical simulations and I will discuss strategies to 
correct for their effects. I also present a brief descrip- 
tion of a third mirelated bias source. 

Keywords binaries: general — binaries: visual — 
methods: numerical — methods: statistical — stars: 
luminosity function, mass function — open clusters and 
associations: individual (NGC 3603, 30 Doradus) 



1 Introduction 

Most initial mass functions (IMFs) of simple stel- 
lar populations are calculated by positioning the ob- 
servational data in a color-magnitude (CMD) or in 
a Hertzsprung-Russell (HRD) diagram. From stellar 
models one extracts the isochrone appropriate for the 
age and metallicity and compares it with the data to 
obtain the evolutionary masses for the individual stars. 
Stars are then grouped into bins with equal width in 
A(logm) and a function (e.g. dn/dm = Ariri^) is fitted 
to the binned data using minimization. 

The above procedure is subject to different biases 
that can influence the derived parameters. In this con- 
tribution, I discuss two of them and present my plans 
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to analyze a third. The first bias is the one induced 
by using bins of equal width when fitting power-laws to 
binned data. Such an effect is more general than its ap- 
plication to the calc ulation of mass functions and wa s 
discussed in Paper I ( Mafz Apellaniz and Ubeda 2005f) . 
The second bias is caused by the existence of unresolved 
multiple systems, either physica l or chance alignments , 
and was discussed in Paper II ( Mafz Apellanid 2008| l. 
The third bias is related to the effect of random uncer- 
tainties in the mass determinations on the computed 
IMF and will be analyzed in a future paper. 



2 Binning biases 

It is well known that when data are binned and a func- 
tion is fitted to the outcome a bias in the derived 
parameters can be (and usually is) present if there 
is a low number of objects in som e of the bins (see, 
e.g. Bevington and RobinsonI 1992t iNousek and Shuel 
19891 ). The problem originates in the strong anticorre- 
lation between the data and the weights in chi-square 
minimization for data w ith P oisson or binomial uncer- 
1995f) . In the case of IMF cal- 



tainties (Wheaton et al. 



culations using equal-width bins biases can be large be- 
cause the values of the slope 7 are typically large and at 
the high end there are usually only a handful of stars. 

There are several ways to deal with this bias. In 
Paper I we introduced the idea of binning the data by 
placing the same number of stars in each bin. Such a 
technique shifts the information on the IMF from the 
number of stars per bin to the positions of the limits 
between bins. Since each bin has then the same weig 



^ If the total number of stars divided by the number of bins is not 
an integer, one has to place one more star in some bins compared 
to others, leading to slightly different weights. Nevertheless, such 
an effect is very small in most cases. 



2 



3.0 



2.5 - 



2.0 



-a 1.5 - 



1.0 



0.5 - 



0.0 - 



★ equal width 

* equal number 



★ * 



10 



30 



100 300 



<A' /bin> 

stars 



Fig. 1 Bias as a function of tiie mean number of stars per 
bin for tiie equal-widtli and tlie equal-number experiments. 



there is no anticorrclation between the data and the 
weights, thus minimizing the biases. 

The idea of using equal- number bins instead of equal- 
width bins was tested in Paper I using Montecarlo sim- 
ulations with different total number of stars (30, 100, 
300, or 1000) and number of bins (3, 5, 10, 30, 50). For 
each combination of numbers of stars and bins 1000 re- 
alizations of a Salpeter IMF with M > 6.3 Mq were 
obtained and the slope of each one of them was mea- 
sured with a x^-fitting algorithm. Then, for each of the 
combinations the bias b was defined as: 



b^< (7fe + 2.35)/afc >, 



(1) 



where k is the realization running index and "fk , <^k are 
the measured IMF slope and uncertainty, respectively, 
for the fcth realization. 

Results are shown in Figure [TJ Fits with equal- 
number bins consistently yield smaller biases than fits 
with equal-width bins. The effect is especially large 
when the mean number of stars per bin is low: for val- 
ues as low as 20, the bias can be larger than 1 when 
equal-width bins are used. On the other hand, even 
when only one star per bin is used, the bias for the 
equal-number case remains smaller than 0.3. Further- 
more, we tested that the uncertainties produced by the 
X^-fitting algorithm indeed correspond to the real un- 
certainties because the standard deviation of the dis- 
tribution of (7fc -f 2.35)/o'/c is very close to 1.0 in all 
cases. 



An interesting side outcome of our Montecarlo exper- 
iments is the determination of the minimum number of 
stars required to measure the slope of an IMF given its 
expected random uncertainty. The lower the number 
of stars, the larger the uncertainty because sampling 
effects become more important. We found the num- 
ber to be ~100 stars for at = 0.2 and ~30 stars for 
CTfe = 0.3. 



3 Unresolved systems 

It is believed that most (if not all) massive stars are 
formed in multiple systems and that only runaway ob- 
jects hav e a relatively low p robability of having no com- 
panions ( Mason et al.1119981 ). Given the typical distance 
to massive young clusters (MYCs) , most of those bina- 
ries will be unresolved even with HST. Also, given the 
distances and core densities of MYCs, chance superpo- 
sitions between cluster members should be a common 
occurrence. An example of the importance of the prob- 
lem is shown in Fig. [21 

In Paper II I modeled the effect of real multiple 
systems and chance alignments on the observed color- 
magnitude diagrams. I also studied how the apparent 
mass function (AMF) differs from the IMF if those ef- 
fects are not taken into account. Below I summarize 
the main results. 



3.1 Real multiple systems 

In the first experiment of Paper II, I modeled the case in 
which all stars in a cluster are unresolved binaries, the 
IMF for individual stars is Kroupa or Kroupa-like, and 
the secondary mass function is flat (hence, for a given 
primary mass, the companion will be, on average, more 
massive than if drawn from an IMF clipped at the pri- 
mary mass). The observed CMD for a Kroupa IMF 
and a comparison between the IMF and the AMF are 
shown in Figs. [3] and [H The observed CMD is widened 
due to the existence of binaries and a population of ap- 
parently ultramassive stars is formed: if the real stellar 
upper mass limit is 120 M0, the apparent one is 182 
M0. The effect on the mass function slope for massive 
stars is small, though, as the AMF and the IMF are 
nearly parallel in Fig. H) 

3.2 Chance alignments of single stars and binaries 

In the second and third experiments of Paper II, I 
treated the effect of chance superpositions due to the 
increased crowding that takes place as the distance or 
central density of a cluster becomes larger. This was 
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done by means of a "fat pixel" approximation that con- 
siders that the Hght from a given star falls on a single 
detector pixel. As crowding increases, each pixel in- 
cludes light from more stars and the observed CMD 
departs more from the original one. The second ex- 
periment analyzed the cases where N — 1, 2, 4 . . . 2048 
single stars are included in each detected source. The 
third experiment did the same using binaries instead 
of single stars. The blue and luminous part of the ob- 
served CMD are shown in Figs. [5] and [S] and the derived 
AMFs are shown in Figs. [7] and [H 

The observed CMD for the second experiment show 
that, as the number of superimposed stars increases, 
the area populated with objects becomes first wider and 
then narrows down. As it happened with the first ex- 
periment, a population of apparently ultramassive stars 
also forms but it takes a relatively large number of su- 
perimposed stars for its number to become significant. 
At the other end, the increase in the number of su- 
perpositions implies that we lose our ability to detect 
low-mass stars. The results for the third experiment are 
similar except that for the case where = 1 (which is 
equivalent to the first experiment) a significant number 
of apparently ultramassive stars already exists. 

With respect to the slope of the AMF for massive 
stars for the second and third experiments, there are 
three regimes. For N small, the change in the slope 
is small (< 0.2), has a weak dependence on N, and is 
easily correctable. For intermediate values of N (~ 30), 
the change in the slope is larger and more dependent on 
A'^ but a correction can still be used. For large values 
of N the situation is hopeless and the true IMF cannot 
be measured from the data. 

3.3 Sample applications: NGC 3603 and R136 

I have applied the results of the experiments above to 
two MYCs in the Local Group, NGC 3603 (Fig. O and 
R136 (Fig.Iini), using HST data (HRC/ACS in the first 
case, PC/WFPC2 in the second one). 

For NGC 3603 chance superpositions are a very small 
effects: most point sources in the HRC image are in- 
deed single stars or unresolved real multiples. This im- 
plies that the massive-star AMF slope (between —1.6 
and —2.0) is indeed very close to the IMF slope. The 
NGC 3603 observations indicate that three stars have 
masses above 120 M©, which is precisely the number 
of apparently ultramassive stars predicted in the first 
experiment if the real stellar upper mass limit is 120 
Mq. Two of t hose three objects are indeed spectro- 
scopic binaries (jMoffat et al.l 120041 Schnurr & Moffat, 
private communication) and there is a relatively separa- 
tion gap where the third one could have an undetected 
(by spectroscopic or visual means) companion. 



For the WFPC2 HST data chance superpositions 
play a more important role: near the core each ap- 
parent point source is expected to be composed of the 
superposition of ^ 2^ individual stars or 2^ bina- 
rie^. Such a large number makes the measurement of 
an IMF slope at the core impossible with the current 
data (the situation could change in the near future with 
new HST data, though). R136 also contains three stars 
with an observed mass of 120 Mq or above. However, 
the existence of those three objects is also compatible 
with them being composed of ^ 100 stars with masse 
below 120 Mq. 

Therefore, the conclusion from the analysis of the 
NGC 3603 and R136 data is that the case for a stellar 
upper mass limit is reinforced and that a value as low 
as 120 Mq for solar metallicity is still possible. 



4 Random uncertainties and the gamma bias 

A third effect that we plan to analyze in the future is 
that of random uncertainties in the determination of the 
IMF. The massive-star IMF should be a monotonically 
decreasing function so one would expect that, given an 
observed mass, the true mass would be more likely to be 
lower than higher than that value because in the parent 
distribution there are more stars with lower than with 
higher masses. Such a "diffusion" from lower to higher 
masses should create an AMF flatter than the real IMF, 
thus creating a "gamma bias" towards higher (lower in 
absolute value) values of the slope, an effect similar to 
the one created by the Malmquist bias for luminosities 
or by the Lutz-Kelker bias for distances. 
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Fig. 2 [top] HST/ACS+WFPC2 3-color mosaic of Pis- 
mis 24, a stellar young cluster located at a distance of 2.5 
kpc. The field size is 2.9 pc x 1.9 pc. [middle] A simulation 
of the same field observed with HST at the distance of the 
LMC. [bottom] A simulation of the same field observed with 
HST at the distance of M33. Note that the brightest appar- 
ent point source in the top frame is actually composed of 
three stars with masses above 50 Mq , of which two of them 
remain unresolved with HST at 2.5 kpc. 



Flat mass-ratio binaries Kroupa 




Fig. 3 [top] Color-magnitude density function for real bi- 
nary systems with a Kroupa IMF and a fiat mass ratio 
shown as a Hess diagram. The black line shows the po- 
sition of the 1 million year isochrone. The function scaling 
is logarithmic and the normalization is arbitrary. 
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Fig. 4 IMF and AMF for real binary systems with a 
Kroupa IMF and a flat mass ratio. The mass is expressed 
in solar units. The IMF is normalized to 1.0 and the AMF 
to 0.5. 
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Fig. 7 IMF and AMFs for chance superpositions of single 
stars with a Kroupa IMF. The mass is expressed in solar 
units. 



Fig. 9 ACS/HRC F550M DRZ frame of NGC 3603. The 
three most massive (WNha) objects are marked and circles 
with radii of 3" and 6" centered on the cluster have been 
drawn. A square-root scale between and 800 counts has 
been used in order to show both the bright and the dim 
stars. The pixel size is 0'.'025, the field size is 29'.'3x25'.'4 
(1172 px X 1016 px), and the vertical direction is 124° East 
of North. 




Fig. 8 Same as Fig.[7]for chance superpositions of binaries. 
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Fig. 10 PC (WFPC2) F555W image of R136. The three 
most massive (WNha) objects are marked and a circle with 
radius of 1" centered on the cluster has been drawn. A 
square-root scale between and 400 counts has been used 
in order to show both the bright and the dim stars. The 
pixel size is 0'.'0455, the field size is 4'.'23x3'.'64 (93 px x 80 
px), and the vertical direction is 37.8° West of North. The 
physical size (1.0 pc x 0.9 pc) is similar to that of Fig. |9l 



Fig. 5 Top left corner of the color-magnitude density functions for A'^ = 1, 2, 4 . . . 2048 chance superpositions of single 
stars and a Kroupa IMF shown as Hess diagrams. The thick black line shows the position of the 1 Ma isochrone and the 
thin lines the evolutionary tracks between and 2 Ma for the initial masses of 25, 40, 60, 85, and 120 Mq. The function 
scaling is logarithmic and the normalization is arbitrary. 
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Fig. 6 Same as Fig. [5] for chance superpositions of binaries. 
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